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Abstract: In arid and semi-arid stream-dominated systems, the temporal variability in groundwater 
recharge has not been widely addressed. Various questions remain about the sources of groundwater 
recharge, its patterns, and the appropriate measuring techniques. Hence, the main objective of the present 
study was to assess the changes that might affect the pattern of groundwater recharge under wetter than 
normal surface water availability. Therefore, the groundwater depth was monitored near a semi-arid 
Mediterranean intermittent stream on the piedmont of the High Atlas Mountains in the mountain 
catchment of the Wadi Rheraya over two hydrological years (2014—2016) with different climate conditions: 
extreme wet and normal conditions. Groundwater recharge was assessed using the episodic master 
recession algorithm. During the two years, the pattern of groundwater recharge was dominated by 
episodic events and by a high seasonality from wet seasons to dry seasons. In the wet year (2014—2015), 
the highest groundwater recharge was recorded following an extreme flood, which deeply replenished 
groundwater. Furthermore, an exceptional steady state of the groundwater depth was induced by a steady 
groundwater recharge rate. For several groundwater recharge events, the assessed recharge had multiple 
sources, mainly from streamflow at the local scale, but possibly from precipitation, underflow, deep 
percolation or irrigation return from the upstream part of the catchment. Local recharge by streamflow 
was likely to be short-lived, and lateral recharge was likely to last longer. Consequently, the episodic master 
recession algorithm estimated the total groundwater recharge that could encompass various sources. In the 
future, more studies and multidisciplinary approaches should be carried out to partition these sources and 
determine their specific contributions. In semi-arid stream-dominated systems, different groundwater 
recharge patterns induced by extreme hydrological events (e.g., wet events) and various potential sources 
of groundwater recharge should be considered when assessing and predicting groundwater recharge. 
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1 Introduction 


In arid and semi-arid basins, numerous plain aquifers are experiencing chronic groundwater 
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depletion because of high groundwater abstraction and/or low groundwater recharge (van Lanen 
and Peters, 2000; Boukhari et al., 2015; Watto and Mugera, 2016; Ashraf et al., 2021). Low 
groundwater recharge is usually due to low precipitation and high evapotranspiration (Flint et al., 
2002; Walvoord et al., 2002). Therefore, in these basins, mountain front groundwater recharge 
linked to streamflow is often considered the primary source of recharge (de Vries and Simmers, 
2002; Manning and Solomon, 2003; Simmers, 2003; Wilson and Guan, 2004; Dahan et al., 2008). 
The infiltration processes of streamflow have been widely analyzed (Ronan et al., 1998; 
Stonestrom and Constantz, 2003; Niswonger et al., 2005; Keery et al., 2007; Rau et al., 2010; 
Lautz, 2012; Clutter and Ferré, 2018; Kurylyk et al., 2019). However, there are few applied 
studies on the processes of groundwater recharge and its temporal variability in arid and semi-arid 
basins (Shanafield and Cook, 2014). 

In general, the behavior of groundwater recharge in arid and semi-arid stream-dominated 
systems is still not well analyzed (Cuthbert et al., 2016). Indeed, when studying groundwater 
recharge in these systems, particularly those with shallow groundwater, various questions might 
need to be addressed. Existing studies have shown that the main source of groundwater recharge 
is from the water losses of nonperennial streams, mainly from floods that bring high amounts of 
water during episodic recharge events (Shentsis and Rosenthal, 2003; Dagés et al., 2008; Cuthbert 
et al., 2019; Shanafield et al., 2021). However, in the case of agricultural practices, substantial 
groundwater recharge could be generated by irrigation returns from diverted streamflow 
(Bouimouass et al., 2020). Therefore, groundwater recharge could be from multiple sources rather 
than a single source. Moreover, the typical pattern of groundwater recharge is characterized by 
episodic recharge, which is related to the intermittency of streamflow and the episodic character 
of floods. In this context, research studies on the variation in the pattern of recharge under 
specific climate conditions are still lacking, e.g., extreme wet events could impact the 
hydrological regime of stream catchments and hence the subsequent groundwater recharge. 
Finally, the water table fluctuation method is one of the most used techniques due to its 
advantages over the other methods of being simple to apply and of providing the most direct 
possible observations and measurements of recharge (Cuthbert et al., 2016; Nimmo and Perkins, 
2018). However, the water table fluctuation method might not always be applicable because of 
some limitations described by Healy and Cook (2002). 

In the present study, a well tapping into a shallow alluvial aquifer near a semi-arid 
Mediterranean intermittent stream was monitored over two hydrological years (2014—2016): one 
with an extremely wet climate (2014—2015) and the other with a typical arid climate (2015-2016). 
Groundwater recharge was assessed over the two years by the water table fluctuation method 
using the episodic master recession algorithm. The main objectives of this study were to (1) 
analyze the changes that may affect groundwater recharge pattern under wetter than normal 
surface water availability; (2) determine the effects of an extreme wet hydrological event on 
groundwater recharge; and (3) explore the adequacy of the episodic master recession algorithm 
for assessing groundwater recharge in semi-arid stream-dominated systems. Thus, this study 
highlighted the complexity of assessing and predicting groundwater recharge in a developed 
semi-arid stream catchment. 


2 Data and methods 


2.1 Study area and data collection 


The study area is located 35 km south of Marrakech City in Morocco, near a semi-arid 
Mediterranean intermittent stream on the piedmont of the High Atlas Mountains in the mountain 
catchment of the Wadi Rheraya (Fig. 1). The area is characterized by a typically Mediterranean 
semi-arid climate. The average annual precipitation (1970-2016) on the piedmont is 350 mm. The 
mountains receive more precipitation (500 mm on average during the period of 2008-2018) both 
as rainfall and snowfall. The mountain catchment of the Wadi Rheraya at the Tahanaout gauge 
station is 225 km? and encompasses the highest peak in North Africa (Jbel Toubkal, with the 
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elevation of 4167 m a.s.l.) (Fig. 1). Streams in the study area are generally intermittent or transient 
types. Part of the streamflow is derived for irrigating olive trees and wheat crops. The alluvial 
unconfined aquifer in the study area comprises alluvial deposits of Neogene and Quaternary age, 
which extend downstream to the plain. The measured transmissivity values vary from 5.0x10> to 
9.0x10 m/s with an average of 6.7x10° m’/s, and the hydraulic conductivity values vary from 
4.6x10~ to 5.5x10 m/s with a mean of 4.4x10~ m/s (Sinan and Razack, 2006). Previous studies 
have presented piezometric maps of the alluvial aquifer showing evidence of groundwater recharge 
from streamflow in the piedmont area (Boukhari et al., 2015; Bouimouass et al., 2020). The 
groundwater is used for drinking purpose and irrigation in dry seasons. 
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Fig. 1 Location of the monitoring well and the mountain catchment of the Wadi Rheraya with its streamflow 
and precipitation gauge stations. 1, Tamatert; 2, Tachedert; 3, Aremd; 4, Matate; 5, Imskerbour; 6, Asni; 7, 
Tahanaout. Groundwater depth is indicated around the monitoring well. 


The present study extended over two years, from July 2014 to August 2016. The measured data 
included daily precipitation in the mountain catchment, streamflow at the Tahanaout gauge 
station, and groundwater depths at a distance of 80 m from the Wadi Rheraya streambed. Daily 
precipitation data over the mountain catchment were measured in the period of 2009-2016 (Fig. 
1; Table 1) in the framework of the Tensift Basin observatory (http://trema.ucam.ac.ma). The 
longest precipitation time series has been acquired on the piedmont of the High Atlas Mountains 
by the Tahanaout gauge station since 1971. 

Daily streamflow has been monitored since 1962 at the Tahanaout gauge station, located in the 
hillslope region, 8 km upstream from the study area (Fig. 1). The main floods were specifically 
measured between 1962 and 2006; afterward, only daily streamflow was measured. 

The groundwater levels were measured at a groundwater monitoring well equipped with two 
automatic measuring probes. A submersible datalogger for real-time water level monitoring 
(TD-Diver probe, Van Essen Instruments, Delft, The Netherlands) based on a pressure sensor 
measured the pressure exerted by both the water column and the atmospheric pressure. Another 
datalogger (Baro-Diver probe, Van Essen Instruments, Delft, The Netherlands) measured the 
real-time atmospheric pressure to account for the changes in the groundwater level induced by the 
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atmospheric pressure alone. The equivalent hydrostatic pressure of the water was used to 
calculate the groundwater level. The groundwater monitoring well was selected to be far from any 
direct pumping influence. First, the monitoring well was in an area surrounded by bare soil and 
natural vegetation. Second, the groundwater level did not show any abrupt drawdowns or 
perturbations that might be induced by pumping. The groundwater depth at the monitoring well 
was 9.00 m at the start of the monitoring (2 July2014), and it varied between a minimum value of 
7.69 m and a maximum value of 10.72 m during the period of 2014—2016. 


Table 1 Precipitation gauges and their measurement periods 


Gauge station Number Longitude Latitude arg Measwrement period 
(m) Start time End time 

Tamatert 1 07°52'48"W 31°08'31"N 1924 18 June 2013 31 December 2016 
Tachedert 2 07°50'56"W 31°09'25"N 2343 1 September 2009 31 December 2016 
Aremd 3 07°55'23"W 31°07'34"N 1940 1 December 2009 31 December 2016 
Matate 4 07°56'28"W 31°09'22"N 1751 19 June 2013 31 December 2016 
Imskerbour 5 07°56'24"W 31°12'25"N 1404 1 September 2009 31 December 2016 
Asni 6 07°58'52"W 31°15'00"N 1170 1 September 2009 31 December 2016 
Tahanaout 7 07°57'47"W 31°17'31"N 1064 1 September 1970 31 December 2016 


The value of the specific yield was estimated to be 0.23 near the study area (Fakir et al., 2021) 
and was the difference between the quasi-saturated water content and the field capacity (Heppner 
et al., 2007). The field capacity value was inferred from water content measurements of the 
alluvial sediment in one hydrological year using soil moisture sensors (Theta Probe Type ML2x, 
Delta-T Devices Ltd., Cambridge, England) at different depths. This estimated specific yield is 
consistent with the value indicated by a previous study of the alluvial aquifer (Sinan, 2000) and is 
comparable to the average value of 0.21 assigned by Healy (2010) to alluvial aquifers. 


2.2 Groundwater recharge assessment using the water table fluctuation method 


In the present study, the water table fluctuation method (Healy and Cook, 2002) was used to 
estimate groundwater recharge by applying the episodic master recession algorithm developed by 
Nimmo and Perkins (2018). The water table fluctuation method stipulates that rises in 
groundwater levels in unconfined aquifers are due to the episodic recharge reaching the water 
table. Groundwater recharge was calculated as the product of two key terms in Equation 1: 


R=AHxS,, (1) 


where R is the groundwater recharge (m); AH is the difference between the peak of the rise of 
groundwater level and the low value of groundwater level estimated for the extrapolated 
precedent recession phase, i.e., the level if no groundwater recharge had occurred; and Sy is the 
specific yield (unitless). The method is adequate for measuring individual and short-term (hours 
or a few days) groundwater recharge events of shallow aquifers with rapid groundwater-level 
fluctuations and for which groundwater drainage can be approximated by graphic extrapolation. It 
can also be applied to assess long-term groundwater recharge rates (seasonal or annual); in that 
case, since groundwater drainage might be significant, what is measured is an estimate of change 
in subsurface storage (m°), referred to as net groundwater recharge (Healy and Cook, 2002). The 
water table fluctuation method considers groundwater recharge in a given period as the sum of the 
episodic recharges calculated from individual groundwater-level rises. 

The main uncertainties of the method could be the frequency of groundwater level 
measurements and the difficulty in determining an accurate specific yield (Crosbie et al., 2005). 
Furthermore, the method excludes situations of high drainage and steady rate of recharge (Healy 
and Cook, 2002). In this study, the water table variation was measured at an hourly time step, 
which is ideal for the water table fluctuation method. The specific yield was estimated near the 
study area (Fakir et al., 2021) as explained in the previous section. 
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The episodic master recession algorithm automatically identifies the main discrete recharge 
episodes and estimates groundwater recharge generated during each of these episodes based on the 
water table fluctuation rates. Three site-specific parameters are required to apply the episodic 
master recession algorithm (Nimmo et al., 2015; Nimmo and Perkins, 2018), which can be inferred 
from the water table data: (1) lag time, the time interval between the occurrence of a water-input 
event and its resulting water table response; (2) fluctuation tolerance, the maximum value of water 
level variation as a function of time, which can be considered system noise rather than a response 
to incoming recharge flux; and (3) master recession curve, a mathematical representation of 
expected water table decline (recession) in the absence of episodic recharge as a function of 
hydraulic head (Crosbie et al., 2005; Heppner and Nimmo, 2005; Nimmo and Perkins, 2018). The 
master recession curve is considered as the relationship between the water table elevation (H) and 
the water table decline rate (dH/dt, where t is the time). 

Based on a first analysis of the observed water table fluctuation, the lag time was equal to 1 d, 
and the fluctuation tolerance was equal to 0.01 m. A lag time of 1 d was justified by the rapid 
reaction of the water table to the hydrological events and the multiday duration of recharge 
episodes. A low value was assigned to the fluctuation tolerance (0.01 m), which was small 
enough to detect all episodes of significant recharge. Both parameters were determined once for 
a given site; consequently, they did not affect episode-to-episode comparisons (Nimmo et al., 
2015). The master recession curve parameters estimated based on the curve fit of dH/dt as a 
function of H, were established over the study period as a linear regression function with a 
coefficient R? of 0.64. Since the groundwater monitoring well was near the wadi and far from 
pumping influence, and as the groundwater table measurements used were free of 
barometric-pressure oscillations, it was considered that the recorded rises in groundwater level 
were due to general recharge (not related to a single source) and the drops in groundwater level 
were attributed to natural discharge. 


3 Results 


3.1 Precipitation analysis 


We performed the seasonal analyses according to the four seasons: autumn (September— 
November), winter (December—February of the following year), spring (March—May) and 
summer (June—August). The annual scale refers to the conventional hydrological year from 
September to August of the following year. Over the mountain catchment, the average annual 
precipitation during the period of 2009-2016 varied from 305 mm at the Asni gauge station to 
595 mm at the Tachedert gauge station (Table 2). The year 2014-2015 was by far the wettest (650 
mm), with a high precipitation gradient across the mountain catchment, while precipitation was 
low (261 mm) and homogeneous in the year 2015-2016 (Fig. 2). 


Table 2 Annual and average precipitation over the mountain catchment of the Wadi Rheraya 


Annual precipitation (mm) 


Gauge station Altitude (m) 
2014-2015 2015-2016 2009-2016 
Tachedert 2343 1093 284 595 
Aremd 1940 860 296 473 
Tamaterte 1924 684 261 - 
Matate 1751 383 198 - 
Imskerbour 1404 600 247 393 
Asni 1170 281 268 305 
Tahanaout 1064 647 275 426 
Average 650 261 - 


Note: -, no data available. 
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Fig. 2 Precipitation isohyetal maps of the mountain catchment of the Wadi Rheraya for the hydrological years 
2014-2015 (a) and 2015-2016 (b) 


The standardized precipitation index (SPI) (McKee et al., 1993) recommended by the World 
Meteorological Organization for meteorological drought monitoring (World Meteorological 
Organization , 2012) has been applied since 1971 using the Tahanaout gauge station data. The SPI 
is calculated as: 

Sp , (2) 
o 
where X; is the precipitation for the station (mm); X is the mean precipitation (mm); and cis the 
standard deviation (mm). The drought severity classes based on the SPI are presented in Table 3. 


Table 3 Drought severity classes based on the standardized precipitation index (SPI) values (McKee et al., 
1993) 


SPI Drought severity class SPI Drought severity class 
>2.00 Extremely wet —1.49 to —1.00 Moderate drought 
1.50-1.99 Very wet —1.99 to -1.50 Severe drought 
1.00-1.49 Moderately wet < -2.00 Extreme drought 
—0.99 to 0.99 Near normal 


For over 46 a of precipitation data in the mountain catchment (Fig. 3), the series was dominated 
by normal precipitation years, constituting 76% of the series, dry years constituting 13% and wet 
years constituting 11%. The year 2014-2015 was classified as the extremely wet year, and the 
year 2015-2016 were categorized as the normal year. 


3.2 Streamflow analysis 


Over the 1962-2016 period, the average annual streamflow of the Wadi Rheraya was 1.55 m?/s. 
The hydrological regime (Fig. 4) shows an important snowmelt contribution in spring and early 
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summer (Hajhouji et al., 2018). The analysis of the main flood events between 1962 and 2006 


showed that these floods were generally rapid (Fig. 5), and 77% of the floods had a base time less 


than 30 h. 
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Fig.4 Average monthly streamflow of the Wadi Rheraya during the period of 1962-2016 
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Fig.5 Variations in the peak flow and time base of the main floods from 1962—1963 to 2005—2006 
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The most frequent flood peak flow rate was between 20.00 and 40.00 m?/s (Fig. 6), but some 
exceptional flow rates could be very high and destructive (e.g., 680.00 m/s in August 1995). 
Normally, floods could occur throughout the year but more frequently in autumn and spring 


(Fig. 6). 
(a) Flood peak flow (b) Season 
3 
Flood peak flow (m*/s) Sanne 
= <20.00 = 20.00-40.00 : , 
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Fig. 6 Frequency of the flood peak flow (a) and seasonal distribution (b) of the main floods from 1962 to 2006 


In the wet year of 2014-2015, the average streamflow was 2.97 m*/s. High-flow events were 
frequent in winter, and the largest flood was recorded on 21 November 2014, with a daily 
streamflow of 48.00 m?/s. The baseflow remained high during 2014-2015 (Fig. 7). In the normal 
year of 2015-2016, the average streamflow decreased to 1.17 m?/s. High-flow events were less 
frequent, and the highest streamflow was 12.93 m/s on 16 February 2016. The baseflow 
decreased dramatically in summer. 
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Fig. 7 Variations in precipitation (a) and groundwater depth and streamflow (b) during the period of 2014-2016 
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3.3 Groundwater depth dynamics 


The groundwater depth dynamics were strongly characterized by fluctuations reflecting the 
alternation of recharge episodes, which coincided with periods of precipitation and floods in the 
Wadi Rheraya, as well as discharge episodes (Fig. 7). The magnitude and duration of these 
fluctuations were very different between the wet year and normal year. 

At the start of the wet year in November 2014, a notable rise in the groundwater depth of 2.53 
m was associated with precipitation and the largest flood on 21 November 2014. The groundwater 
depth remained high for several months. In spring, the groundwater depth rose again following a 
flood at the end of March 2015 and remained at a stable level for almost 1.5 months until the end 
of May 2015. From June to September 2015, groundwater experienced a long and smooth 
recession. In the normal year of 2015-2016, the groundwater depth fluctuations occurred more 
rapidly than those in the previous year. After the rise of the groundwater depth in spring, the level 
of the groundwater depth decreased significantly in summer when the streamflow in the Wadi 
Rheraya became very low. 


3.4 Groundwater recharge variation 


From July 2014 to August 2016, 13 groundwater recharge episode events (E1 to E13) were 
identified by the episodic master recession algorithm (Fig. 8; Table 4), eight in the wet year and 
five in the normal year. 

In the wet year of 2014-2015, the highest groundwater recharge (0.66 m) occurred in 
November 2014 (E2; Fig. 8), following an exceptional hydrological event characterized by the 
highest recorded flood (48.00 m?/s) and precipitation (334 mm). The sixth recharge episode (E6; 
Fig. 8), which lasted 44 d, corresponded to the period when the groundwater depth was stable, 
which could be explained by a steady rate of groundwater recharge. In this case, the water table 
fluctuation method developed for episodic events could not accurately calculate groundwater 
recharge. However, using the same master recession curve as the other events, groundwater 
recharge of this episode was estimated to be 0.06 m (Table 4). 
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Fig. 8 Results of the episodic master recession analysis showing the groundwater recharge episodes (E1—E13) 
during the period of 2014-2016. MRC, master recession curve. 


The cumulative estimated episodic groundwater recharge of the wet year was 1.19 m (63% of 
the total groundwater recharge for both years). The winter recorded 70% of the annual 
groundwater recharge. It was mainly generated by E2 on 22 November 2014 (up to 84%). 
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Groundwater recharge occurred over 245 d, or 68% of the wet year. 

In the normal year of 2015-2016, the cumulative episodic groundwater recharge was 0.67 m 
(37% of the total groundwater recharge for both years). Up to 72% of groundwater recharge 
occurred in the wet year, with most occurring in spring (46%). The hydrological events in summer 
did not generate groundwater recharge, leading to a decrease in the groundwater depth. 
Groundwater recharge processes lasted only 131 day, or 38% of the normal year. 


Table 4 Characteristics of the groundwater recharge episodes inferred from the episodic master recession during 
the period of 2014—2016 


Groundwater ae _ Cumulative Total Number of Peak 

Year recharge Start time End time mae peared groundwater precipitation precipitation streamflow 
episode recharge (m) (mm) days (d) (m°/s) 
El 17 Sep 2014 10 Oct 2014 24 0.05 0.05 23.14 3 9.87 
E2 24 Oct 2014 24 Dec 2014 62 0.66 0.71 334.19 19 48.02 
E3 4 Jan 2015 1 Feb 2015 29 0.08 0.79 52.19 8 8.46 
= E4 11 Feb 2015 26 Feb 2015 16 0.06 0.85 13.90 6 2.30 
z E5 8 Mar 2015 6 Apr 2015 30 0.14 0.99 89.62 10 9.41 
° E6 14 Apr 2015 27 May 2015 44 0.06 1.05 50.73 7 4.17 
E7 9 Jul 2015 26 Jul 2015 18 0.02 1.07 30.98 7 2.85 
E8 6 Aug 2015 27 Aug 2015 22 0.06 1.13 36.72 8 2.58 
E9 20 Sep 2015 9 Nov 2015 5l 0.19 1.32 55.43 10 3.73 
Ka E10 7 Jan 2016 14 Jan 2016 8 0.03 1.35 1.68 1 1.88 
g Ell 9 Feb 2016 9 Mar 2016 30 0.14 1.49 55.21 10 12.93 
R E12 16 Mar 2016 1 Apr 2016 17 0.14 1.63 46.36 6 3.05 
E13 27 Apr 2016 21 May 2016 25 0.17 1.80 19.07 4 2.48 


For the whole study period, the cumulative groundwater recharge was 1.80 m. The total 
duration of groundwater recharge was estimated to be 375 d, which was 48% of the observation 
period. The duration of groundwater recharge events was longer than the duration of the 
corresponding precipitation events (Table 4). This may indicate that groundwater recharge was 
mainly due to streamflow. In addition, there was a low correlation between groundwater recharge 
and peak streamflow (r=0.19, R?=0.03) or average streamflow (r=0.30, R?=0.03) of each event; 
these coefficients were calculated without considering the exceptional second event (E2) because 
of its large difference from the others. 


4 Discussion 


The study area is located near a semi-arid Mediterranean intermittent stream that is generally 
characterized by a high snowmelt contribution from the High Atlas Mountains, which are 
responsible for higher streamflow and more floods in spring and early summer. In the mountain 
catchment of the Wadi Rheraya, the snow regime is a main characteristic, as shown by several 
previous studies (e.g., Hanich et al., 2022). The annual contribution of snow to streamflow was 
assessed from 15% to 50%, depending on the meteorological conditions of the year (Boudhar et 
al., 2009; Marchane et al., 2013). According to the SPI, the climate on the piedmont is 
characterized by 13% of the years being considered dry and 11% being considered the wet years, 
and the remaining 76% of years are considered normal. 

The present study focuses on the changes that may affect groundwater recharge mechanisms in 
the wet year. In both the wet year (2014—2015) and the normal year (2015-2016), the majority of 
groundwater recharge occurred during the wet seasons. In the wet year, the very rainy winter, with 
its associated streamflow and floods, was the highest groundwater recharge period (almost 70% 
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of the total), while in the normal year, as winter had lower precipitation, the spring season with its 
floods issued from snowmelt was the highest groundwater recharge period (46%). Therefore, the 
time of greatest groundwater recharge was during the periods with higher precipitation, higher 
streamflow or more floods when the amount of water could saturate the underlying sediment and 
generate high amounts of recharge (McCallum et al., 2014). In summer, there was little (in the 
wet year) or no recharge (in the normal year). Seasonality is a main characteristic of the temporal 
distribution of groundwater recharge under a semi-arid climate. Seasonality was explained by the 
high moisture of the sediment in the wet season (Fakir et al., 2021), which generated higher 
hydraulic conductivities and fewer unsaturated pores than that in the dry season characterized by 
dry sediment due to water scarcity. 

In the wet year, two additional features were observed. First, the largest flood on 22 November 
2014 was classified as an extreme event with a return period estimated to be approximately 20 a 
(Zkhiri et al., 2017). The amount of surface water generated by this event led to the above normal 
groundwater level rise and subsequently to a very high recharge amount. The latter represented 
60% of the recharge recorded in the wet year. The alluvial aquifer was significantly replenished. 
Second, in spring 2015, an exceptional steady state of groundwater depth occurred during almost 
1.5 months due to a steady rate of groundwater recharge. In this case, the rate of groundwater 
recharge was constant and equal to the rate of groundwater drainage. This groundwater recharge 
regime seemed unusual in semi-arid mountain catchments where the typical pattern of 
groundwater is mainly characterized by episodic recharge. For that regime, groundwater recharge 
was not episodic, and its estimation fell beyond the scope of the water table fluctuation method 
(Healy and Cook, 2002). However, it was possible to use the episodic master recession algorithm 
to estimate this recharge. Other techniques could be used in such conditions, but they require 
further investigations and aquifer characteristics (Eaton, 2019). In semi-arid mountain 
catchments, such sustained groundwater recharge can only be linked to exceptionally high and 
constant streamflow. 

The recorded duration of groundwater recharge events, often several weeks, was longer than 
that of the corresponding precipitation events. It was also longer than the duration of floods that 
are generally rapid (77% of floods were less than 30 h). Furthermore, a low correlation existed 
among recorded groundwater recharge, precipitation and streamflow. Therefore, groundwater 
recharge near the streambed seemed to be due not only to local recharge from streamflow but also 
to upstream lateral recharge from several potential sources, including underflow, deep percolation 
from the upper reach of the catchment or irrigation return from upstream. Local recharge by 
streamflow is likely to be short-lived, and lateral recharge at a slower rate is likely to last longer. 
Therefore, the episodic master recession algorithm estimated the amounts of total recharge that 
probably integrated various components. Hence, the time lag used in the episodic master 
recession algorithm should not be considered a representative characteristic of a unique source of 
recharge. 

This fact emphasizes the importance of combining other methods (e.g., hydrochemistry, 
isotopes, etc.) to assess each component of the total recharge. Consequently, in semi-arid 
stream-dominated systems, the variation in groundwater recharge is rather erratic because it is 
controlled by numerous parameters. Among them are the abovementioned multiple sources of 
recharge and their time scale, sediment moisture conditions (Batlle-Aguilar and Cook, 2012; 
Schwartz, 2016; Fakir et al., 2021), soil depth above the water table (Cao et al., 2016), sediment 
permeability (Shanafield et al., 2012; Villeneuve et al., 2015), climate variability (Moeck et al., 
2020) and land use changes (Lasagna et al., 2020). Therefore, assessing and predicting the 
variability of groundwater recharge in such developed semi-arid mountain catchment is complex 
and needs to consider the various parameters controlling the processes of groundwater recharge. 
In addition, long-term and large-scale monitoring could better characterize the temporal 
variability and the catchment-scale heterogeneity of groundwater recharge. 
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5 Conclusions 


Water table fluctuations were monitored near a semi-arid Mediterranean intermittent stream in a 
semi-arid mountain catchment over two hydrological years (2014—2016) with contrasting climate 
conditions: a wet year and a normal year. In general, groundwater recharge was episodic and 
presented a high seasonality from wet seasons to dry seasons. Extremely wet periods brought a 
large amount of water that introduced an unusual groundwater dynamic and hence recharge rates. 
An extreme hydrological event induced a very large rise in the groundwater depth. Additionally, 
an unusual steady state regime of the groundwater depth was observed, induced by a steady 
groundwater recharge rate. These conditions strongly replenished the alluvial aquifer. 

Groundwater recharge near the semi-arid Mediterranean intermittent stream on the piedmont of 
the High Atlas Mountains was due to local recharge mainly from streamflow and probably from 
heavy precipitation of the wet year. However, the duration of some groundwater recharge events 
was longer than the recorded floods or precipitation events, showing the importance of lateral 
recharge from various potential sources, probably combining underflows or irrigation returns 
from a large part of the stream-dominated systems. Therefore, the episodic master recession 
algorithm allowed us to estimate the total groundwater recharge that could encompass various 
sources. More research should be carried out at large temporal and spatial scales to segregate 
these sources and determine their respective contributions. This is vital for assessing and 
predicting groundwater recharge in such a developed semi-arid mountain catchment, especially 
considering the threat of climate change. 

In arid and semi-arid regions, wet conditions are rare, and their return period is generally long. 
Consequently, despite their importance for groundwater recharge and aquifer replenishment, they 
cannot counterbalance the effects of droughts and mismanagement that have severely affected 
groundwater resources. 
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